General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



/ v 7 7 yv 



DtPART M ENT Of MATH EMATICS 
UNIVERSITY Ol HOUSTON HOUSTON. TEXAS 


(NASA-CF-14774C) USEES GUIDE UHHIE/FTCC N76-24927 

VEFSICN; FROGPAf* ECCU P E N I A I ICN S1GFXT 
EOC-LAPSYS VERSION CF 7 H H L E (Houston Univ.) 

19 p HC $3.50 CSCL 09E Unclas 

G 3/6 1 4 2113 





•Rfo 'V R 1 0 
l arih nr v rva r i 

IJNuJ 

CONTRACT "luS 


FOR 

111 . 

£R 

- >-12777 



HOUSTON. TEXAS 77004 



USERS GUIDE UHMLE/RTCC VERSION 


March 1976 


William A. Coberly 
University of Tulsa 


Consultant to 


Department of Mathematics 
University of Houston 


Report if 52 


NAS-9-12777 



USERS GUIDE UHMLE/RTCC VERSION 


♦ Step 1 Load 

• Step 2 EXEC 
* 

Medium 

CARD 

CARD 

TAPE 

DISKFILE 


DISKFILE 


Program tape (MLETAP « 013856). 

PGM * 0LDTAP 

INPUT 

Variable Description Format 

•NCHS No. of channels 

KCH(I) ,1=1, . . . >NCHS Channel list (1712) 

KFILE CK/RT tape file no. (15) 

CK/RT on file KFILE of Fortran unit 1. 

OUTPUT 

Statistics from CK/RT tape on Fortran unit 12 
Mean vector 

Class 1 

Covariance (symmetric ) 

• 

« 

Mean vector 

Class NCLASS 

Covariance 

NCLASS - No. of classes (16) 

. m 

fortran unit 13 


• Step 3 EXEC PGM = MLF. 


(UHMLE) 



INPUT 


CARD 

RUNID 


(10X,A4) 

DISK 

M 

No. of classes (unit 13) 

(16) 

CARD 

N 

No. of channels 

(10X,I/0) 

CARD 

L 

buffer size 

II 

CARD 

FMT(18) 

format of temp, data set 

(18A4) 

CARD 

KFILE 

file no. of image data 

(15) 

CARD 

SAMKEY 

(=1) 

(10X,I10) 

CARD 

ISTART 

start line 

absolute 

II 

CARD 

IST0P 

stop line 


CARD 

ISKT.P 

line skip (~ ,1 no skip) 

It 

CARD 

JSTART 

start pixel 

relative 

II 

CARD 

JST0P 

stop pixel 


CARD 

JSKIP 

pixel skip (= 1 no skip) 

II 

CARD 

NCH0UT 

number of channels from tape 

VI 

CARD 

NCHLST( ) 

channel list 

(10X.16I2) 

CARD 

NFLDS 

If NFLDS = 0 

no. of fields input 

skip field def cords and process 
all of the test site 

(lOX.IlO) 

CARD 

FID , NV , MINLIN , HAXLIN 

(A8 > 2X, 315) 

'card 

IF (NV+1) 

line coord. 

field 1* 

(1115) 

CARD 

JF (NV+1) 

pixel coord. 

(1115) 

CARD 

FID, NV, MINLIN, 

,MAXLIN 


CARD 

IF (NV+1) 

field NFLDS 

it 

CARD 

JF (NV+1) 




* first field should be definition of test site all fields 
'$$$$$$$$’ are deleted from the temporary data set. 


with FID 


TAPE Image data on file KFILE Universal 

Unit 1 

CARD YUNIT Temp stat unit (= 12) 15 

DISK Stats on YUNIT (8F10.5) 

CARD* IFLAGA 
IFLAGM 
IFLAGS 
M0DES 
ITLIM 

T0LA (One Card) 

T0LM 

T0LS 

EPSA 

EPSM 

EPSS 

ITRPNT (5I5,6F5.3,I5) 

* This card may be repeat as many times as necessary 

OUTPUT 

DISK Temporary image data set unit 23 FMT 

DISK MLE output STATS on unit 11 (8F10.5) 

• Step 4 EXEC NEWTAP 

INPUT 

CARD NCHS ,KCH(NCHS) (1712) 

CARD KFILE CK/RT tape file no. (15) 

TAPE' Old CK/RT tape on unit X 

DISK Stats from MLE unit 11 (8F10.5) 

OUTPUT 

New CK/RT tape on unit 2 


TAPE 
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•SIGEXT Processor 


The SIGEXT processor is a version of the UH1ILE program, developed at the 
University of Houston, which has been Integrated into the Univac Exec II LARSYS 
system. The program uses initial signatures from a training site and unlabeled 
data from a recognition site and computes maximum likelihood estimates of the 
mixture distributions of the unlabeled sample. That is, estimates are found 
for the mixture proportions the mean vectors y^, and the covariance 

matrices E^, for i * 1, . ..,m. For a theoretical discussion of the algorithms 
see [1,2,3.]* 

Two options are available: 

(1) General form : There is no constraint on the movement of the subclass 

signatures in the Iteration process. That is, there is no a priori trans- 
formation formaly assumed. 

(2) Affine form : It is assumed that there exists a diagonal matrix A and 
a vector b defining an affine transformation Ax + b» Maximum likelihood 
estimates for A and b are found and the output signatures are computed as 
follows : 


* Ay ± + b, = A T E^A, i ■= 1,2,..., m 

For a more detailed engineering description and program documentation see [4] . 

The following preliminary user’s guide is intended to conform to the for- 
mat of 15]. Hopefully, any unresolved terms, formats, or references will be 


found there. 
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Input files 

An MSS data tape In Universal or LARSYS format must be assigned to logical unit CL 
If initial statistics are input from tape, logical unit A must be assigned 
to the tape. If the statistics are »put from a Module Stat Deck, then a Fastran 
file must be assigned to logical unit A, 

2. Output files 

The output statistics may be saved on file 2 of logical unit A or punched 
on cards. 

3. Scratch files 

A drum file is used as an intermediate image data file. 

A temporary Fastran file, logical unit \ is required when estimating the 
ALPHA parameter vector only, 

4. Card Input 

All required input card-types are described below. The actual formats for 
card input are described in Section .3 of thfe LARSYS Users Guide (see ref [5]). 

4.1 PROCESSOR CARD 

The SIGEXT processor must be activated by the f $SIGEXT f processor card. 


Keyword * 


Function 

$ SIGEXT 


Directs the LARSYS system 
executive monitor to load 
the signature extension pro- 



cessor and begins input of pro 



cessor control cards. 


■k 

Briefly, the KEYWORD must begin in Column 1 and the parameter field 
begins in Column 11 and ends in Column 72. 



4.2 SPECIAL SYSTEM DECKS 


The training statistics may be input by means of the 'Module STAT deck*. 
The LARSYS system deck formats are describe in SECTION 3 of the LARSYS Users 

Guide • , 

4.3 PROCESSOR CONTROL CARDS 


The following control cards are the complete set recognized by the signature 
extension processor. All options available to the user are exercised by use 
of the appropriate processor control card. If a default condition is specified 
the control card is optional, with the processor using the default parameters 
if the control, card is not input. If no default condition is specified, the 


control card is mandatory input to the processor. Ordering the sequence of 


processor control cards is unnecessary, with the exception of the '*END*' card 


and ’SEND*- card . The ? *END*' card must follow the last processor control 
card the ' $END* * card must follow the last field definition card defining an 


area to be processed. 

KEYWORD 

MODULE 


PARAMETER (S) 


(DEFAULT : training 
subclass statistics 
will be read from the 
input file, SAVTAP) 


FUNCTION 

Indicates to the processor 
that the training subclass 
statistics will be card- 
input. The 'MODULE STAT 
Deck' must immediately 
follow this 'MODULE' 
control card. See Section 
3.14.1 ref I5J for further 
detail on this deck. 



KEYWORD 

PARAMETER (S) 

FUNCTION 

OPTION 

STATS 

Training subclass statistics 


(DEFAULT: no 

will be printed out (mean vector 


training sub- 

and covariance matrix) for each 


class statistics 

subclass. 


printout. 


0PTI0N 

PUNCH 

Output statistics from the 


(DEFAULT: no stat 

SIGEXT Processor will be 


deck is punched 

punched on a ’M0DULE STAT 

DECK*. 

SUBCLASS 

K 1' K 2 K m 

K i a ’ ;s inte 8 ers 

The subclass numbers listed 
will be the set of subclasses 


where 1SK.560, 
t 

used by the processor. 


m * no. of sub- 



classes in train- 

This set must be a 


ing statistics. 

subset of the training sub- 


(DEFAULT: All 

classes on the SAVTAP file- 


the training sub- 

K 1 * K 2 * • * ‘ ,K m are Integer 


classes will be 

numbers designating sub- 


used.) 

classes as they occur on 



the statistics file. 

HED1 

Any 60 alphanumeric 

The Processor will use these 


^ characters beginning 

60 characters as the first line 


in col, 11. 

of the page heading for every 


(DEFAULT: standard 

page of printout produced. 


heading, "LYNDON B. 

(Blanks are Included as 


JOHNSON SPACE 

’characters’) 


CENTER’) 




KEYWORD 

ALPHA 


CHANNEL 


PARAMETER (S) 

ApAg, . . . , 

or 

n*a 1 ,k*a 2 ,a 3 

where N and K 

are arbitrary re- 

pltltion factors and 

A^ are decimal numbers *. 

nt Is the number of 

subclasses. The A 's 

1 

are automatically normalized 
so that Z A^ “ 1 0 

(DEFAULT: If no ALPHA card 

is included then A . * — for 

1 m 

each 1 — l^.o* ,m») 

TRAINING * N. , . . . ,N , TEST = L. , . * . ,L 
1 n’ 1* n 

Channels N. ..... N are 

1’ * n 

selected from the input 

training statistics file 

and channels L* , . . . ,L are 
1 n 

selected from the input image 
data file, 

(DEFAULT: all channels are 
selected from both stat file 
and the image data file. Warning 
these files must have the same 
number of channels) . 


FUNCTIONS 
These entries are 
the initial values 
of the mixture pro- 
portions in the MLE 
procedure 



KEYWORDS 

PARAMETER (S) 

FUNCTION 

HED2 

Any 60 alpha-* 

The processor will 


numeric charac- 

use these 60 characters 


ters beginning 

as the second line oj; 


in col. 11. 

the page heading, for 


(DEFAULT: stan- 

every page of print-out 


dard second line, 

produced (blanks are 


’HOUSTON, TEXAS') 

included as characters') 

COMMENT 

Any 60 alpha- 

The processor will use 


numeric charac- 

these 60 characters 


ters beginning in 

as the third line of 


column 11. (DEFAULT: 

the page heading, for 


all blanks in third 

every page of print-out 


line of page heading) 

produced (Blanks are 
included as 'characters') 

DATE 

MM DD YY any (up to) 12 

The processor places 


alphanumeric characters 

these 12 characters 


beginning in column 11. 

in the upper right- 


(DEFAULT: the current 

hand corner of each 


date, at the time of the 

page of print-out. 


run, will be placed in 

Normally used to 


upper right-hand corner 
of each page of print- 
out) . 

specify a calendar date. 

*END* 


This control card term- 


(DEFAULT: none) 

inates the input of pro- 
cessor control cards, and 
initiates the input of 
data card(s) defining 
area(s) to be processed. 
The field definition 
card(s) for the area(s) 
to be processed must imm- 
ediately follow the ,;k ENE 


control card. 



KEYWORD 


FUNCTION 


$END* -none- This control card terminates 

the input of field definition 
card(s) defining area(s) to 
be processed. This control 
card must immediately follow 
the last field definition 
card. 

4.4 FIELD DEFINITIONS 

Area(s) to be processed are communicated to the signature extension pro- 
cessor by using the 'field definition 1 data card, described in Section 3.1.3, 

of Ref. [5], These cards contain.the scan line and sample coordinates for the 

vertices of tb.o polygon 3haped area over which signature extension is* performed. 

At least one 'field definition' card must be in the run deck, immediately following 

the '*EKD*' control card. As many 'field definition' cards as there are area(s) may 

be input. The areas specified on the 'field definition' card(s) must be available 

on the MSS data file that is input to the processor. 

4.5 Iteration Control Cards 

KEYWORD PARAMETER 

0PC0DE 1 or 2 

(DEFAULT 1) 


FUNCTION 

If 1 Is chosen, the 
general MLE procedure 
is usedo If 2 is 
chosen, then the Ax + b 
form is used. 



KEYWORD 


IF LG A 


IFLGM 


IFLGS 


IFLGAM 


IFLGBM 


M0DES 


PARAMETER 


(DEFAULT Ij « 1 
for j * 1, . . . ,m) 


(DEFAULT = 1 
for j ** 1,. . . *n) 

X l»' * * *\ 

(DEFAULT I. * 1 
J 

for j * 1, ... fa ) 


I 

(DEFAULT 1=1) 

I 

(DEFAULT I * 1) 

1 or 2 
(DEFAULT 1) 


FUNCTION 

If I. “0, then the mixture 
proportions for the j 
class is fixed for this 
iteration phase. If 

I. * 1, then the mixture 

J th 

proportions for the j 

class is updated. 

If I . * 0, then the mean 

J th 

vector for the j class 

is fixed for this iteration 

phase. Otherwise, the j C ^ 

mean vector is updated. 

If 1^ * 0, then the j*'* 1 
covariance matrix is fixed 
for this Iteration phase. 
Otherwise the j covariance 
matrix Is updated. 

If I = 0 then the diagonal 
matrix A is fixed for this 
phase of the Ax + b iteration 
Otherwise the matrix A is 
updated. 

If I = 0,then the vector b is 
fixed for this phase of the 
Ax + b iteration. Otherwise, 
the vector b is updated. 

If 1 is chosen, then the 
full covariance matrix is 
used. Otherwise, only a 
diagonal covariance matrix 


EPSA 

0 < X <; 2 
real number 
(DEFAULT X = 1.0) 

X Is the step size para- 
meter for the proportion 
iterations. 

EPSM 

0 < X £ 2 
real number 
(DEFAULT X - 1.0) 

X is the step size para- 
meter for the mean vector 
iterations 

EPSS 

0 < X £ 2 
real number 
(DEFAULT X * 1.0) 

X is the step size < ti- 
me ter for the covar.i ; vice 
matrix iteration. 

EPSAM 

0 < X £ 2 
real number 
(DEFAULT X * 1-0) 

X is the step size para- 
meter for the A matrix 
iteration. 

EPSBM 

0 < X S 2 
real number 
(DEFAULT X » 1.0) 

X is the step size para- 
meter for the b vector 
iteration. 

T0LA 

X > 0 

real number 
(DEFAULT X * .001) 

Iteration tolerance on pro- 
portions. The maximum change 
in proportions from one 
iteration to the next over 
all classes is compared to 
T0LA in the iteration termi- 
nation decision. 

T0LM 

X > 0 

real number 
(DEFAULT X = .5) 

Iteration tolerance on mean 
vectors. The maximum change 
in mean vector components over 
all channels over all classes 
from one iteration to the next 
is compared to T0LM in the 
iteration termination decisions. 

T0LS 

X > 0 

real number 
(DEFAULT X * 1. ) 

Iteration tolerance on covarianc 
matrices. The maximum change in 

covariance matrix elements over 
all channels over all classes 
from one iteration to the next 
is compared to T0LS in the iter- 
ation termination decision* 


T0LAM 


T0LBM 


ITRPNT 


ITLIM 


*END* 


$END* 


X > 0 

real number 
(DEFAULT X = .05) 


X > 0 

real number 
(DEFAULT X * .05) 


I 

integer 

(DEFAULT I * OX 


I 

integer 

(DEFAULT I =» 10) 


none- 


-none- 


Iteration tolerance on 
the matrix A in the Ax + b 
option. The maximum change 
in the elements of the diagonal 
matrix A from one iteration to 
the next is compared to T0LAM 
in the iteration termination 
decision. 

Iteration tolerance on the 
vector b in the Ax +b 
option. The maximum change in 
the elements of the vector 
from one iteration to the next 
is compared to T0LBM in the 
iteration termination decision. 

Iteration print control. 

If I * 0 only the results of 
the first and last iterations 
for each iteration phase are 
printed. Otherwise, all 
iterations are printed. 

Limit on number of iterations 
permitted. The iteration pro- 
cess is terminated for a given 
iteration phase after ITLIM 
iterations have occurred re- 
gardless of tolerance comparison 

This- card follows each set of 
iteration control parameters 
which define an iteration phase. 

Control is returned to the 
monitor from the SIGEXT pro- 


cessor. 



5. Deck Setup for LARSAA (SIGEXT) 


there may be any number of 
iteration control blocks. 





*END* 


Iteration Control Cards (Phase 1) 


$END* 


field definition cards 


*END* 

Module STAT DECK (Optional) 
Processor Control Cards 


$SIGEXT 


$ XQT 
TRW Z 

IN Z 
TRW Z 


LARSAA 


program 

Loads program LARSAA from tape 


@ XQT CUR 


0 ASG I (temporary file) 


@ ASG C = Image data tape 


@ ASG A * SAVTAP 


@ ASG Z = Program tape 


. If Unit A is input tape then 
File 1 is input stats 
File 2 may be used to save output 
stats. 

. If stats are created by $STAT,$ISOCLS 
or read from a MODULE STAT DECK then a 
position file should be assigned. 


@ Z RUN 




6. Flow of program 


SIGEXT^> — _JJo ^ Branch to another processor 


Enter SIGEXT.SETUPA routine to 
read processor control parameters 
and set defaults 


read data indicated by field definition ' 
cards from tape to disk file 

i 

reads initial stats for Indicated channels 
& classes from SAVTAP 
~ 1 

Print initial st ats 

V_ 

Enter iteration supervisor UHM SUP 



Enter SETUPB routine to read iteration 
control parameters 


Perform computations required by the 
iteration control parameters 


return to SIGEXT routine write output 
stats if required punch output stats 
if required 


return to LARSAA Monitor 









